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Abstract 

The recently developed structure model that uses the generator coordinate method to perform 
configuration mixing of angular-momentum projected wave functions, generated by constrained 
self-consistent relativistic mean-field calculations for triaxial shapes (3DAMP+GCM), is applied in 
a systematic study of ground states and low-energy collective states in the even-even magnesium 
isotopes ^'^"^"Mg. Results obtained using a relativistic point-coupling nucleon-nucleon effective 
interaction in the particle-hole channel, and a density-independent (5-interaction in the pairing 
channel, are compared to data and with previous axial IDAMP+GCM calculations, both with a 
relativistic density functional and the non-relativistic Gogny force. The effects of the inclusion of 
triaxial degrees of freedom on the low-energy spectra and E2 transitions of magnesium isotopes 
are examined. 

PACS numbers: 21.10.Ky, 21.10.Rc, 21.30.Fe, 21.60.Jz, 27.30. +t 
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I. INTRODUCTION 

In the first part of this work [l| the simple mean-field (single-reference) implementation 
of the framework of relativistic energy density functionals (EDF) has been extended to in- 
clude long-range correlations related to restoration of symmetries broken by the static mean 
field and to fluctuations of collective coordinates around the mean-field minimum. A model 
has been developed that uses the generator coordinate method (GCM) to perform con- 
figuration mixing of three-dimensional angular-momentum projected (3DAMP) relativistic 
mean-field wave functions, generated by constrained self-consistent calculations for triaxial 
nuclear shapes. This calculational framework can be used to perform detailed studies of low- 
energy collective excitation spectra and corresponding electromagnetic transition rates. The 
particular implementation of the relativistic 3DAMP+GCM model has been tested in the 
calculation of spectroscopic properties of low-spin states in ^"^Mg, in comparison with data 
and with the results of the recent work of Ref. [2I, where a similar 3DAMP+GCM model 
has been developed, but based on non-relativistic Skyrme triaxial mean-field states that are 
projected both on particle number and angular momentum, and mixed by the GCM. We 
note that, very recently, a new 3DAMP+GCM model with particle number projection has 
been implemented, based on the non-relativistic Gogny force 

In this work we apply the relativistic 3DAMP+GCM model to a systematic study of 
ground states and low-energy collective states in the even-even magnesium isotopes 20-40Mg. 
The low-energy structure of magnesium nuclei has attracted considerable interest in the last 
decade, both experimental and theoretical. In particular, the sequence of isotopes ^°~^°Mg 
encompasses three spherical magic shell numbers: = 8, 20, 28, and, therefore, presents an 
excellent case for studies of the evolution of shell structure with neutron number, weakening 



of spherical 
inversion" 



shell closures, disappearance of magic numbers, and the occurrence of "islands of 



J]. Following the pioneering measurement of the transition rate B{E2; Of 2f) 
in the neutron-rich nucleus ^^Mg ^| that confirmed a large deformation of this nucleus 
indicated by the low excitation energy of the 2f state jo], extensive experimental studies of 
the low-energy structure of Mg isotopes have been carried out at RIKEN 0,3, MSU IJQ, 
GANIL 13 1 and CERN QQ- 

In addition to numerous theoretical studies based on large-scale shell-model calculations 



16 



2l|, the self-consistent mean- field framework, including the non-relativistic Hartree- 
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Fock-Bogolibov (HFB) model with Skyrme 22| and Gogny forces 23|, and the relativistic 



mean- field (RMF) model 2J, |25|, as well as the macroscopic-microscopic model based on a 



modified Nilsson potential 26|, have been used to analyze the ground-state properties (bind- 
ing energies, charge radii and deformations) and low-lying excitation spectra of magnesium 
isotopes. Of course, to calculate excitation spectra and electromagnetic transition rates it 
is necessary to go beyond the mean-field approximation and include dynamic correlations 
related to the restoration of broken symmetries and to fiuctuations of collective coordinates. 
Based on the IDAMP+GCM (axial symmetry) framework, studies of low-energy spectra of 



specific Mg isotopes have been performed using non-relativistic models with Skyrme 
and Gogny 231] forces, as well as relativistic density functionals 
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3Q|. 



In Section [ITl we present a brief outline of the relativistic 3DAMP+GCM model used 
in the present analysis. Section IIIII describes a study of low-lying collective states of the 
even-even magnesium isotopes ^°~^''Mg. A brief summary and an outlook for future studies 
are included in Section I 



II. THE 3DAMP-KGCM MODEL 

The IDAMP-I-GCM calculational framework, restricted to axially symmetric nuclei, has 
recently been extended to include triaxial shapes. 3DAMP+GCM models have been devel- 
oped, based on the self-consistent Hartree-Fock-Bogoliubov approach with Skyrme forces [2I 
and the Gogny force 3|. Starting from relativistic energy density functionals, we have imple- 
mented a model for configuration mixing of three-dimensional angular-momentum projected 
(3DAMP) relativistic mean-field wave functions, generated by constrained self-consistent cal- 



nuclear shapes. The details of the model and the numerical tests are 



I31I 1. Here we only outline the basic features of the model that will be 



culations for triaxia 
described in Refs. [l 
used in the study of low-lying states in even-even magnesium isotopes 20-40Mg. 

In the 3DAMP+GCM framework the trial angular-momentum projected GCM collective 
wave function |\E';^*'^), an eigenfunction of J^, Jz with eigenvalue J(J + l)/z^ and Mh, reads 



where a = 1,2, ■■ ■ labels collective eigenstates for a given angular momentum J, and q is 
the generic notation for the deformation parameters /3 and 7. The projection of the angular 



momentum J along the intrinsic z-axis K takes only non-negative even values, and Pl/m 
denotes the angular-momentum projection operator: 

2 7+1/" 

p'mk = J dnDi:An)R{n) . (2) 

Q denotes the set of three Euler angles: {0, 6*,^/;}, and dQ = d(j) sin 6 dO dtp. D'l,j^{VL) is the 
Wigner D-function, and the rotational operator reads R{VL) = ^'^'t'J^Q^sJvQ^i'J^ ^ The set of 
deformed intrinsic wave functions \^{q)) is generated by imposing constraints on the axial 
q2o and triaxial 522 niass quadrupole moments in self-consistent RMF-I-BCS calculations. 

The weight functions fi^{q) in the collective wave function Eq. ([I]) are obtained from 
the solution of the Hill- Wheeler-Griffin (HWG) integral equation: 

I dq' J2 [-^KK'iQ, q') - Ei^A'iQ, q')] = 0, (3) 

K'>0 

where and ^ are the angular- momentum projected GCM kernel matrices of the Hamil- 
tonian and the norm, respectively jlj] . The solution of Eq. determines both the energies 
and the amplitudes fa^{q) of collective states I^E';^*^) with good angular momentum. 
The center-of-mass correction to the total energy of the state is calculated in the zeroth 
order of the Kamlah approximation. 

Since the weight functions fi^{q) are not orthogonal and cannot be interpreted as col- 
lective wave functions for the deformation variables, the collective wave functions ga{i) are 
determined from the eigenstates of norm overlap kernel 

9i^) = Y.9t4{i)- (4) 

k 

These functions are orthonormal and 

i 

where the sum is over i = {K, q}. The coefficients (7^° are solutions of the following equation 

Y,nigf^ = Eigi\ (6) 
I 

which is equivalent to Eq.(l3]). The matrix "H^^ is determined by the angular-momentum 
projected GCM kernel matrix of the Hamiltonian 



where and wl are the non-vanishing eigenvalues and eigenvectors of the norm overlap 
kernel respectively. 

The B{E2) value for a transition from an initial state (Jj,aj) to a final state {Jf,af) is 
calculated from 

2 



B{E2; Ji,ai Jf,af) 



2Ji + 1 

where the reduced matrix element is defined by 



^{'Jf^<if\\Q2\\Ji,qi) 



{Jf,qf\\Q2U,q^) = (2J/ + 1) 5^ 



fiK' 



J f 2 Ji 
-Kf /i K 



(8) 



(9) 



with f^^{q) = (—1)"^/^ ^{q) fo^^ K < 0. The spectroscopic quadrupole moment for the 
state (J^) is defined by the expression 



IGtt 



J -J 



J 2 J 



x(2j + i)(-i)^+^'5^ I " ^ " \{HqJ)\Q2,P^>>Kmq^))■ (10) 

The matrix elements of the charge quadrupole operator Q2fj, = €.J2p''^p^2fj.{^p) are calcu- 
lated in the full configuration space. There is no need for effective charges, and e simply 
corresponds to the bare value of the proton charge. 



III. LOW-LYING STATES IN MAGNESIUM ISOTOPES: RESULTS AND DIS- 
CUSSION 



As in the first part of this work we use the relativistic point-coupling interaction PC- 



Fl 



32| in the particle-hole channel, and the corresponding density-independent 5-force in the 



particle-particle channel. The parameters of the PC-Fl functional and the pairing strength 
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constants Vn and Vp have been adjusted simultaneously to the nuclear matter equation of 
state, and to ground-state observables (binding energies, charge and diffraction radii, surface 
thickness and pairing gaps) of spherical nuclei 32j, with pairing correlations treated in the 
BCS approximation. In particular, the pairing strength parameters for neutrons and protons 
are K = -308 MeV fm^ and Vp = -321 MeV fm^, respectively. 

Parity, D2 symmetry, and time-reversal invariance are imposed in the constrained mean- 
field calculation of the binding energy map of an, in general triaxial, even-even nucleus. To 
solve the Dirac equation for triaxially deformed potentials, the single-nucleon spinors are 
expanded in the basis of eigenfunctions of a three-dimensional harmonic oscillator (HO) in 
Cartesian coordinates, with Ng^ = 8 major shells for ^°~^^Mg and A^sh = 10 for ^^~^°Mg. 
These numbers of oscillator shells are sufficient to obtain a reasonably converged mean-field 
potential energy surface [3l|. The HO basis is chosen isotropic, that is the oscillator 
parameters = hy = = = ^Jh/muJQ in order to keep the basis closed under rota- 



tions 
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341 . The oscillator frequency is given by fkOQ = AlA . The Gaussian-Legendre 



quadrature is used for integrals over the Euler angles 0, 9, ip in the calculation of the norm 
and hamiltonian kernels. The choice of the number of mesh points for the Euler angles in 
the interval [0, tt] is: = = 8, and Ng = 12. In the 3DAMP-FGCM calculations of 
^'^Mg it has been shown that, because of very few level crossings as function of deformation, 
redundancies appear very quickly in the norm kernel when more states are added to the 
nonorthogonal basis |l|, ^ . The generator coordinates are, therefore, chosen in the intervals 
< /3 < 1.2 and < 7 < 60°, with steps A/3 = 0.2 and A7 = 20°, respectively. Moreover, 
eigenstates of the norm overlap kernel with very small eigenvalues ri'^, / n"^^^ < C are removed 
from the GCM basis. With the cutoff parameter ( = 5x 10"^ for ^°"^^Mg, and C = 1 x 10~^ 
for ^^"^"Mg, fully converged results are obtained for all low- lying states with J < 6. 

In Fig.[T]we plot the self-consistent RMF-I-BCS mean-field, and the corresponding angular 
momentum projected ( J'^ = 0''"), energy curves (PEC) for the even-even magnesium isotopes 
2o-40]y[g^ as functions of the axial deformation /3 (7 = 0). One might notice an interesting 
evolution of the mean- field PECs from a spherical shape at magic neutron number = 8, 
through pronounced prolate shapes, coexistence of oblate and prolate shapes, and again 
to a spherical shape at = 20. Increasing further the neutron number from A^ = 20 
to A^ = 28, the mean-field minima become markedly prolate. Furthermore, the effect of 
angular momentum projection can be inferred from a comparison with the corresponding 
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{J'^ = 0^) PECs in the right panel of Fig. [T] In particular, in the neighborhood of the 
spherical minimum the J = PECs of ^°''^^Mg are very soft with respect to f3. In other 
isotopes the deformed minima become deeper after projection. 

Figure [2] displays the total ground-state dynamical correlation energies of Mg isotopes, as 
a function of the number of neutrons. As shown in the figure, -Ecorr consists of a rotational 
energy correction AEj=o that results from the restoration of rotational symmetry 

AEj=o = Ej^oiM - EuFiPm) , (11) 

and the correlation energy gained by GCM configuration mixing 

AEgcm = ^(0^) - ^j=o(/3o) . (12) 

and Po denote the axial deformation parameters at the minima of the mean-field and 
the (J'^ = 0"*") angular-momentum projected PECs, respectively (cf. Fig. [1]). Ecott shows a 
strong dependence on shape and shell structure. It is large for deformed mid-shell nuclei, 
with a maximum of ~ 4 MeV at = 14, and is drastically reduced (~ 1 MeV) for the two 
isotopes with the neutron magic numbers N = 8 and = 20. Projection on angular momen- 
tum J = , that is, the rotational energy correction AEj=q constitutes the dominant part of 
the total dynamical correlation energy. This is generally valid for a great majority of nuclei, 



as it has been shown in the global study of quadrupole correlation effects [35|, performed 
with GCM configuration mixing of axially symmetric Skyrme-Hartree-Fock+BCS states, 
with the two-point topological Gaussian overlap approximation for angular-momentum pro- 



jection. As also shown in Ref. 35|, in Fig. [2] one notices that the correlation energy Ai^QCM 
gained from configuration mixing of different deformed states is of the order of several hun- 
dreds keV, and not very sensitive to nuclear shape and shell structure. AE'gcm is, in fact, 
composed of two parts: a potential term that is negative and in size comparable to the 
correlation energy induced by angular momentum projection, and a kinetic part (energy of 
the zero-point vibrational motion) that is positive and cancels to a large extent the potential 
term 36|. 

The excitation energies of the states 2f and in ^*^~^°Mg, calculated using the 
IDAMP+GCM model with the relativistic density functional PC-Fl, are compared in Fig. [3] 
to the available data and the prediction of the IDAMP-I-GCM calculation based on the non- 
relativistic HFB framework with the Gogny force Both models yield excitation energies 
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of the 2f and 4^*" states in reasonable agreement with data and, on the average, the values 
obtained with PC-Fl are 10 — 30% lower than those calculated with the Gogny interaction 
DIS (except for ^^Mg). This is due to relatively weak neutron pairing correlations in the 
present calculation, that lead to an increase of the corresponding moment of inertia for the 
yrast states. As noted in our previous study of ^^Mg in Ref. the excitation energies of 
yrast states increase when the pairing strength parameters Vn/p are adjusted to the pairing 
gaps determined from empirical odd-even mass differences in this particular mass region. 
Both calculations preserve the N = 8 magic number and with PC-Fl also at = 20 a 
pronounced shell closure is obtained, whereas the model based on the Gogny force predicts 
a much lower excitation energy of the 2f state in ^^Mg, in better agreement with data. One 
might notice, however, that both models predict the state in this nucleus at energies far 
above the experimental value. The N = 28 shell closure disappears in both calculations, 
and ^°Mg is predicted to be prolate deformed. 

The corresponding B{E2;0t 2^) (e^fm^) values in 20-40j^g 

are shown in Fig. |H 

IDAMP+GCM calculations, both the present one using the functional PC-Fl and that 



based on the Gogny force [23[, yield results in reasonable agreement with data except, of 
course, PC-Fl at and in the neighborhood of the neutron number A^ = 20. Since the Gogny 
force DIS predicts an axially deformed ground state for '^^Mg, the corresponding B(E2) 
value for the transition Of — 2^ is much closer to the experimental value, compared to the 
calculation with PC-Fl which yields a spherical ground state at A^ = 20. The functional 
PC-Fl, together with the density-independent 5- force {Vn = —308 MeV fm^ and Vp = —32 



MeV fm"^) predicts indeed a very small B{E2) value for this transition in "^^Mg. In Ref. 31 1 



it has been suggested that a better adjustment of pairing strength parameters and eventu- 
ally the inclusion of triaxiality, that is the 7 degree of freedom, could improve the results 
for ^^Mg. Already in the IDAMP-I-GCM axial calculations we have verified that, by adjust- 
ing the pairing strengths specifically to the empirical pairing gaps around ^^Mg (five-point 
formula): Vn = —465 MeV fm^ and Vp = —350 MeV fm^, the calculated transition rate 
increases to B{E2;0i — ?■ 2^) = 313.5 e^fm^. To have a consistent model, however, in the 
remaining calculations of this work we will continue using the original pairing strengths that 
were adjusted simultaneously with the parameters of the PC-Fl effective interaction in the 



particle-hole channel 



32|. 



A measure of collectivity of the lowest excited states in magnesium isotopes can be 
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obtained by comparing the experimental and calculated B(E2; Oi — )■ 2^) with the prediction 
of an empirical formula based on the liquid-drop model (LDM) 38| . 

B{E2 : Oi 2+)sys. = 6A7Z^A-^-'^^E-\2t) . (13) 

This comparison is shown in Fig. [5l In the upper panel the B(E2) values calculated with 
the LDM formula are compared to data, whereas in the two lower panels they are compared 
to the results of the IDAMP+GCM calculations with the functional PC-Fl and with the 
Gogny force DIS. The excitation energies E{2f) (in MeV) that appear in the LDM expression 
Eq. ( IT3|) . correspond to the experimental values and those calculated with PC-Fl and Gogny 
DIS, respectively. One notices a very good agreement between data and the B(E2) values 
predicted by the LDM formula. Based on the recently measured E(2i) values for ^°Mg: 
1598(10) keV, and ^^Mg: 660(6) keV, Eq. ([HD predicts the corresponding B{E2;0f 2|) 
values 368.9(23) e^fm^ and 595.4(54) eHm^, respectively. The IDAMP+GCM calculation 
based on the PC-Fl functional yields somewhat smaller B(E2) values for the Oi — ?> 2^ 
transition in ^OMg (332 e^fm^) and ^^Mg (460 e^fm^). 

In Fig. [6] we plot the spectroscopic quadrupole moments of the states 2^ and 4^ in 
2o-40]y[g^ calculated using the IDAMP+GCM model with the relativistic density functional 
PC-Fl, and compared to the corresponding values based on the non- relativistic HFB frame- 



work with the Gogny force 23|] . One might notice a very good agreement between the results 
of the two model calculations, with the exception of '^°Mg. In the lower panel the calculated 
ratios Q^^'^^{Af)/Q^^'^^(2'l) are compared to the value that corresponds to a rigid axial rotor 
with A' = 0, that is ~ 1.27 . In ^^Mg both models predict a very small value of Q^^'^'^{2^), 
and this gives rise to an exceptionally high value of Q^^'^'^{4:i)/Q^^'^'^{2f) that does not fit 
the scale of the vertical axis. This result indicates that there is a large contribution from 
nonzero- -fT components in the yrast band of ^^Mg. Large deviations from the axial rotor 
value are also predicted for ^°Mg and "^^Mg. For the isotopes 22,24,28,32-40]Yjg i^q^^]^ models 
yield Q^^'^^{4i)/Q^^^^{2f) quite close to that of rigid axial rotor. Note that this is also 
true in "^^Mg, for which the calculation based on the Gogny force yields a deformed ground 
state, whereas this state is spherical in the present axially symmetric calculation using the 
functional PC-Fl. In both calculations, however, the states 2^ and 4^ are prolate deformed. 

To examine the infiuence of triaxiality, that is, of including the 7 degree of freedom on 
the spectroscopic properties of low-lying states in magnesium isotopes, we have performed 
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full 3DAMP+GCM calculations using the relativistic functional PC-Fl. In Figs. [7| and [8] we 
display the resulting self-consistent RMF+BCS triaxial quadrupole binding energy maps of 
the even-even ^''~^°Mg isotopes in the /3 — 7 plane (0 < 7 < 60°), and the corresponding 
angular-momentum J'^ = 0^ projected energy surfaces. All energies are normalized with re- 
spect to the binding energy of the absolute minimum, the contours join points on the surface 
with the same energy (in MeV). In general the inclusion of the triaxial deformation degree 
of freedom reduces considerably the barriers separating axially prolate and oblate minima 
in the well-deformed isotopes 22,24,34-40]Yjg ^^jg^ notice that the angular-momentum 

J"^ = 0+ projected energy surfaces of ^^~^^Mg are rather soft both in /3 and 7. 

The low-energy excitation spectra and collective wave functions are calculated as solutions 
of the Hill- Wheller- Griffin integral equation for each angular momentum, and thus take into 
account fluctuations of the collective coordinates /3 and 7 around the mean-field minima. 
For the sequence of isotopes ^°~^°Mg Figs. |9] and [10] display the probability distributions 
\g;l\'^ of the collective wave functions Eq. (j4]) in the /3 — 7 plane, for the states Oi and 
2f (both the K = and K = 2 components). It appears that 20,30,32^ 

are spherical in 

the ground state, whereas all the other isotopes are prolate deformed and the ground-state 
deformation is especially pronounced in heavier Mg nuclei. The first excited state 2^ is 
prolate deformed in all Mg nuclei, even in ^^Mg. In several isotopes, most notably in ^^Mg 
and ^°Mg, the collective wave function of the state 2f contains sizeable admixtures of the 
K = 2 component. This can be seen more clearly in Fig. [TT] where, after integrating the 
probability distributions over /3 and 7, in the upper panel we plot the relative weight of the 
K = component in the collective wave functions of the 2f states of magnesium isotopes 
2o-40]y[g_ rjj^g softness toward triaxial shapes is especially pronounced in ^^Mg, ^^Mg, and 
^°Mg. The contribution of the K = 2 component in the wave functions of 2^ will generally 
affect the calculated B(E2) values for transitions to the ground state. In the lower panel 
of Fig. [11] we show the differences between the B{E2;2i — t- 0^*^) values calculated in the 
full 3DAMP-FGCM and the axial IDAMP-t-GCM models, normalized to the ID values. A 
marked effect of i^'-mixing is found not only in ^^Mg, but also in some heavier isotopes 
including ^^Mg. 

Finally, a quantitative comparison between the axial IDAMP + GCM and the full 
3DAMP + GCM calculations for ^o-^OMg, based on the relativistic functional PC-Fl, is 
presented in Table [T] The ground-state energies Egg (in MeV), excitation energies of the 2f 
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and 4^ states (in MeV), and B{E2 J,; J — t- J — 2) values (in e^fm^) for the lowest states with 
J = 2"*", 4"'' in magnesium isotopes are included in the table. In general the inclusion of the 
7 degree of freedom leads to the lowering of the binding energies of low-lying states and to 
an increase of the calculated B(E2) values. The latter is particularly prominent in ^^Mg, in 
which the 3DAMP + GCM yields an enhancement of ^ 25% for the B{E2; 0+ 2+). Es- 
pecially interesting is the case of "^^Mg, which shows a pronounced lowering of the excitation 
energies of 2f and 4^, whereas the binding energy of the ground-state, being spherical, is not 
influenced by the inclusion of triaxial shapes. These excitation energies are, however, still 
far above the experimental energies and even though the B{E2; 0^ — )■ 2^) value increases 
by ~ 10%, it is about a factor three smaller than the empirical value. However, when the 
pairing strength parameters are adjusted specifically to the empirical pairing gaps around 
^^Mg (five-point formula): Vn = —465 MeV fm^ and Vp = —350 MeV fm^, the calculated 
transition rate increases to B{E2; Of — t- 2f) = 330.1 e^fm^, in rather good agreement with 
data. This results shows the importance of a more detailed study of pairing correlations in 
^ 20 neutron-rich nuclei. 



IV. SUMMARY 

The very successful framework of relativistic energy density functionals has mostly been 
used on the mean-field level to describe ground-state properties of medium-heavy and heavy 
nuclei. When considering applications, however, it is important to develop EDF-based 
structure models that go beyond the static mean-field approximation. Detailed predictions of 
excitation spectra and transition rates necessitate the inclusion of correlations related to the 
restoration of broken symmetries and to fluctuations of collective variables. In recent years 
several new models have been developed that extend the relativistic EDF-based approach 
and perform the restoration of symmetries broken by the static mean field and take into 
account fluctuations around the mean-field minimum. This is relatively simple in the case 
of axial symmetry, that is, when only one collective coordinate is considered Q , but 
such models become much more involved, technically complicated, and computationally 
demanding when possible triaxial shapes are taken into account. 

In Refs. 31] and [l| we have implemented and tested a new model that uses the gen- 



erator coordinate method (GCM) to perform configuration mixing of three-dimensional 
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angular-momentum projected (3DAMP) relativistic mean-field wave functions, generated 
by constrained self-consistent calculations for triaxial nuclear shapes. In the present study 
this calculational framework has been used to analyze the influence of triaxiality on the 
low-energy collective excitation spectra and the corresponding electric quadrupole transi- 
tion rates of even-even magnesium isotopes ^°~^°Mg. The self-consistent solutions of the 
constrained RMF-I-BCS equations have been obtained using the relativistic point-coupling 



interaction PC-Fl 32] in the particle-hole channel, and a density- independent (5-force in the 
particle-particle channel. Since the low-energy spectra of ^°~^''Mg were previously investi- 
gated in the axial IDAMP-I-GCM model based based on the non-relativistic HFB framework 



with the Gogny force 23|, in the first instance we have performed axial ID calculations and 



compared the results with data and those obtained in Ref. [23[. In general, a good agreement 
has been obtained between the results of the two model calculations, except for ^"'■^^Mg. The 
low excitation energy of 2^ and the large B{E2; Of 2f) indicate that the neutron rich 
nucleus ^^Mg is deformed, even though the number of neutrons equal the "spherical magic 
number" = 20. The data are reproduced reasonably well by the IDAMP+GCM model 
based on the Gogny force, which yields a deformed ground state for ^^Mg. The present axial 
calculation, on the other hand, predicts a spherical /3-soft ground state for '^^Mg, although 
the lowest excited states 2f and 4f are calculated to be prolate deformed. The correspond- 
ing B{E2;0f — 7- 2^) is much smaller than the experimental value. Both models predict 
prolate ground states for heavier Mg isotopes, including the A^ = 28 nucleus ^°Mg. 

To analyze the effect of triaxiality and K-mixing on the low-energy structure of Mg 
isotopes, we have also performed a full 3DAMP+GCM calculation based on the relativistic 
density functional PC-Fl and a density-independent 6 pairing interaction. When compared 
with the IDAMP+GCM results, it is noted that the inclusion of the 7 degree of freedom 
leads to the lowering of the binding energies of low-lying states and to an increase of the 
calculated B(E2) values in deformed isotopes. In several isotopes a pronounced degree of 7 
softness and K-mixing is predicted for the yrast states. The effect is strongest in ^^Mg, in 
which the 3DAMP + GCM yields an enhancement of ^ 25% for the B{E2; Of 2^). Even 
in the triaxial case the functional PC-Fl preserves the spherical shell closure at A^ = 20, i.e., 
it predicts a spherical ground state for ^^Mg. The excitation energies of the states 2f and 
4:f in this nucleus are lowered considerably with respect to the axial case, but they are still 
much higher than the experimental values. Correspondingly, the calculated B{E2; Of — t- 2^) 
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is about a factor three smaller than the empirical value. It is noted, however, that when the 
pairing strength parameters are adjusted specifically to the empirical pairing gaps around 
^^Mg, the calculated transition rate increases to B{E2; 0^ 2f) = 330.1 e^fm"^, much closer 
to the available data. 

In future studies the 3DAMP + GCM model based on relativistic density functionals 
will be applied to the description of shape transitions and shape coexistence phenomena in 
medium-heavy and heavy nuclei. We also plan to compare the results of full 3D angular- 
momentum projection and GCM configuration mixing, with those obtained in the recently 
developed model for the solution of the eigenvalue problem of a five-dimensional collective 
Hamiltonian for quadrupole vibrational and rotational degrees of freedom, with parame- 
ters determined by constrained self-consistent relativistic mean-field calculations for triaxial 
shapes [39 1. 
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TABLE I: The ground-state energy Egg (in MeV), excitation energies of the 2f and 4^^ states 
(in MeV), and B{E2 ■\; J - 2 J) values (in e^fm'') for the lowest states with J = 2+,4+ in 
magnesium isotopes. Results obtained in the axial IDAMP+GCM calculation are compared with 
those of the full 3DAMP+GCM model. 
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FIG. 1: (Color online) Self-consistent RMF+BCS mean-field (left panel), and angular-momentum 
projected 0''" potential energy curves (PEC) (right panel) of even-even magnesium isotopes, as 
functions of the axial deformation parameter (3. To plot all the curves in the same figure, the 
PECs of 20-28]y[g j^^^^g ^^^^ shifted by -75, -55, -35, -20, and -10 MeV, respectively. The position 
of the minimum of each PEC is indicated by a red dot. 
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FIG. 2: (Color online) Total ground-state dynamical correlation energies Eqo„ of Mg isotopes, as 
a function of the number of neutrons. -Ecorr is the sum of the rotational energy correction IS.Ej=q 
and the energy gained by configuration mixing AE'qcm- 



19 



FT 



> 5 

X 3 
2 

X 

liJ 1 


7 

> 6 
S 5 
X 4 
V 3 

LU ^ 
1 





□ 

O 



■ 8 8 

J i L 



1 — ' — \ — ' — I — ' — 
■ Exp. 
□ HFB(Gogny) 
O RMF(PC-FI) 



O Q I 



■BBBo 

J I I i L 




8 12 16 20 24 28 
Neutron number 



FIG. 3: (Color online) Excitation energies of the states 2^ and 4^ in ^'^"^''Mg, calculated using the 
IDAMP+GCM model with the relativistic density functional PC-Fl, are compared to available 
data U, Q, S^l and the results of the IDAMP+GCM calculation based on the non-relativistic 



HFB framework with the Gogny force 
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FIG. 4: (Color online) B{E2;0f 2f) (e^fm^) values in '^^~^^Mg, calculated using the 
IDAMP+GCM model with the relativistic density functional PC-Fl, are compared to available 



data 



3, 0, y, Q, Q, Q 



and the results of the IDAMP+GCM calculation based on the non- 



relativistic HFB framework with the Gogny force 
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FIG. 5: (Color online) B[E2) values calculated with the LDM formula Eq. (I13p for the transition 



OJ" — )■ in Mg isotopes, are compared to data 



flfl 
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371] in panel (a), and to the results 



of the IDAMP+GCM calculations with the functional PC-Fl in (b), and with the Gogny force 
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FIG. 6: (Color online) Spectroscopic quadrupole moments of the states 2f and 4^ in '^^ ^"^Mg, 
calculated using the IDAMP+GCM model with the relativistic density functional PC-Fl, and the 



corresponding values based on the non-relativistic HFB framework with the Gogny force [23|] (upper 
panel). The calculated ratios Q^^^'^{4f)/Q^^^'^{2f) are compared to the value that corresponds to 
a rigid axial rotor with K = (lower panel). 
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FIG. 7: (Color online) Self-consistent RMF+BCS triaxial quadrupole binding energy maps of the 
even-even ^''"^^Mg isotopes in the /3 — 7 plane (0 < 7 < 60*^) (left panel), and the corresponding 
angular- momentum J'^ = 0^ projected energy surfaces (right panel). All energies are normalized 
with respect to the binding energy of the absolute minimum, the contours join points on the surface 
with the same energy (in MeV). 
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FIG. 8: (Color online) Same as described in the caption to Fig. [7] but for the isotopes '^''Mg. 
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FIG. 9: (Color online) Probability distributions I^^P of the collective wave functions Eq. ^ in the 
/3 — 7 plane, for the the states of Oj*" and 2^ (both the K = and K = 2 components) of ^"^"^^Mg. 
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FIG. 11: (Color online) Upper panel: relative weight of the K = component in the collective 
wave functions of the 2f states of magnesium isotopes ^'^~^*^Mg. Lower panel: differences between 
the B{E2;2f 0|) values calculated in the 3DAMP+GCM and the IDAMP+GCM models, 
normalized to the ID values. 
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